Modification of the Turbulent Energy Cascade by Polymer Additives 
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By tracking small particles in the bulk of an intensely turbulent flow, we show that even a very small con- 
centration of long-chain polymers disrupts the usual turbulent energy cascade. The polymers affect scales much 
larger than their physical size, from the dissipation range to the inertial range. The effect depends strongly on 
the polymer concentration. While the dissipative-scale statistics change continuously as the polymer concen- 
tration is increased, the inertial-range energy transfer rate is only altered by the polymer additives when the 
concentration is above a threshold (approximately 5 parts per million by weight for the polymer we used). 



PACS numbers: 47.27.Jv,47.27.Gs,47.57.Ng,47.50.-d 

Minute amounts of long-chain flexible polymers added to 
a fluid can strongly modify flow properties. In a turbulent 
wall-bounded flow, for example, they lead to the extraordi- 
nary phenomenon of drag reduction 111 |2[ |j, |4j, |5|] . These ef- 
fects may be qualitatively attributed to the stretching of poly- 
mer molecules by the velocity gradients in the flow. Although 
progress has recently been made in understanding drag re- 
duction at a wall, comparatively little is known about the 
action of polymers far from the boundaries of a turbulent 

flow JS Umm. 

Fluid turbulence is inherently far from equilibrium: energy 
continuously passes into and out of the system. In a New- 
tonian fluid like water, where molecular viscosity provides 
the sole mechanism for the dissipation of energy, the clas- 
sical Richardson-Kolmogorov cascade hypothesis [10] states 
that energy is injected into the flow at large length and time 
scales, transferred to smaller and smaller scales without loss, 
and finally dissipated at the smallest scales where viscosity 
acts. The rates of energy injection ej, energy transfer ex, and 
energy dissipation en, are therefore equivalent in Newtonian 
fluid turbulence. Long-chain polymer molecules, which tend 
to coil up like balls of thread in their equilibrium state, can 
be stretched by the straining of the fluid flow to many times 
their equilibrium length and store elastic energy. In a turbu- 
lent flow, the polymer molecules will stretch and recoil in the 
fluctuating flow field, and this process will dissipate kinetic 
energy due to interactions between the monomers of one poly- 
mer molecule and between the polymers and the fluid. The 
addition of polymers to a turbulent flow therefore provides a 
new route by which kinetic energy can flow out of the tur- 
bulent phase. The previously unaddressed but fundamental 
question is then the relationship between e/, er, and en for 
turbulence in polymer solutions. 

In this Letter, we show that a very small concentration of 
long-chain polymer molecules strongly modifies the turbulent 
cascade of energy from large to small scales. Our results in- 
dicate that the effect of polymers on the energy cascade can 
be divided into two regimes depending on polymer concen- 



tration. For small concentrations, the energy transfer rate er 
is unchanged and only the viscous dissipation rate erj is re- 
duced due to the additional dissipation mechanism provided 
by the polymers. For large concentrations, however, we find 
that not only e o but also the apparent energy transfer rate er 
decreases. The observations cannot be explained by current 
theories. 

Newtonian turbulence is described by a single nondimen- 
sional parameter, the Reynolds number, which compares the 
strength of inertial driving forces to viscous damping forces. 
A large Reynolds number implies that a wide range of length 
and time scales participate in the energy cascade, where the 
statistical properties of the flow are expected to be universal 
and independent of the driving mechanism. This universal 
regime is known as the inertial range. A polymer solution 
is additionally characterized by both the polymer concentra- 
tion and the Weissenberg number Wi, which compares r p , the 
relaxation time of a single polymer, to the fastest flow time 
scale. In turbulent flows, it is defined as Wi = r p /T n , where 
t v is the Kolmogorov time scale. When the Weissenberg num- 
ber is less than a critical value (of order unity), the polymer 
molecules are generally in their coiled state and will be pas- 
sively advected by the flow. When the Weissenberg number 
is larger than this value, the polymers will be stretched by the 
flow and may modify it . 

Earlier experimental investigations of the interactions of 
bulk turbulence and polymers have generall y ei ther forced 
the turbulence through a boundary layer ifTH Il2l [l3i H3| . or 
have been performed at relatively low Reynolds number {(^fj^, 
where the turbulence was not fully developed and it was dif- 
ficult to quantify the effect of polymers on turbulence, par- 
ticularly for inertial-range quantities. In a water flow be- 
tween counter-rotating disks with raised vanes, it was ob- 
served that the energy injection at the disks remained con- 
stant when adding polymers into the flow 01511 . A previous 
experiment with polymers in our apparatus, however, showed 
that the acceleration statistics of the flow were strongly af- 
fected Jsl ISD . Numerical simulations of isotropic turbulence 
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with polymers have also been performed at small Reynolds 
number using model equations but they are 

very difficult due both to the nature of the equations and to 
numerical instabilities 12011 . 

In this Letter, we report results from Lagrangian parti- 
cle tracking experiments conducted in a water flow between 
counter-rotating disks with vanes. Flow properties were mea- 
sured by tracking [ 21] the simultaneous motion of hundreds of 
nearly neutrally buoyant 33 /mi fluorescent polystyrene tracer 
particles, excited by a high-power pulsed Nd:YAG laser de- 
livering up to 90 W and recorded with three Phantom v7.1 
CMOS cameras from Vision Research, Inc. The polymer used 
was an 18 x 10 6 a.m.u. molecular weight polyacrylamide 
(Poly sciences 18522) with an equilibrium radius of gyration 
of 0.5 /xm, a fully stretched length of 77 /im, and a relaxation 
time of t p — 43 ms [9;]. While ou r appar atus allows us to 
reach very high Reynolds numbers 11221 l23l |24 12511 . the size 
and flexibility of the polymer molecules makes them prone 
to tearing in very intense turbulence; we therefore only con- 
sider Reynolds numbers where our results are not affected by 
polymer degradation fl E3[. In our experiments, the Weis- 
senberg and Reynolds numbers are coupled: based on the 
smallest turbulent time scale, the Weissenberg number ranges 
from Wi = 1.2 to 6.0. The Taylor-microscale Reynolds num- 
ber R\ = yJlhu'L/v ranges from 200 to 350, where u' is 
the root-mean-square turbulent velocity, L is the largest length 
scale of the turbulence, and v is the kinematic viscosity. We 
varied the polymer concentration from (pure water) to 20 
parts per million by weight (ppm). Note that, in order to make 
a quantitative comparison with Newtonian fluid turbulence, 
the Reynolds numbers we report for experiments with poly- 
mer solutions are those measured in pure water before poly- 
mers were added to the flow. 

To measure the modification of the turbulent energy cas- 
cade by the polymer additives, we use the second-order trans- 
verse Eulerian structure function Djvjv(?") = ((8 r u) 2 ), which 
measures the difference in velocity over a separation r, as 
a probe of the scale-by-scale properties of the cascade. For 
the transverse structure function, velocities are measured or- 
thogonal to r. The analogous longitudinal structure function 
Dll(t) is defined with the velocities taken along r. We show 
here only measurements of Z?Arjv(r); our results are equiva- 
lent for DLL(r). D jv\/v (r) has three distinct scaling regimes. 
In the small-scale dissipation range (r <C rj), for isotropic tur- 
bulence, 
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At large scales (r ^ L), DjvAr(r) saturates at twice the veloc- 
ity variance. At intermediate scales, in the so-called inertial 
range (77 <C r <C L), the classical Kolmogorov theory ifioll 
predicts that 
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where C2 = 2.13 ± 0.22 is a well-known universal constant 
determined from previous experiments [26]. We note again 



(b) 



FIG. 1: (color online) The effect of polymer concentration on the 
Eulerian structure functions. The Reynolds number of the water flow 
(before adding polymers) is R\ = 350 (corresponding to Wi = 6.0 
with polymers) and the Kolmogorov length scale is r) w = 84/im. 
The structure functions compensated by (a) the inertial-range scal- 
ing prediction and (b) the dissipation-range scaling prediction. The 
dashed line in (b) indicates et(0) as measured in the water flow. 



that in Newtonian turbulence, the energy transfer rate ex and 
the energy dissipation rate eo are the same. 

We show in Fig. Q] the measured DNN{r) for different 
polymer concentrations at R\ — 350, where R\ was mea- 
sured without polymers and the corresponding Kolmogorov 
length scale is r\ w — 84/im. In Fig. Ufa), we plot DNN(r) 
compensated by the Kolmogorov scaling prediction (i.e., 
[(3/4)D A r A r(r)/C 2 ] 3/2 /(e T (0)r), where e T (0) is the energy 
transfer rate measured from the water data), so that a plateau 
at unity indicates inertial-range scaling. We observe that at 
small concentrations (<j> < 5 ppm in our experiments), the 
shape of DNN(r) remains the same: the curves return to 
the same plateau value in the inertial range, while the ex- 
tent of the apparent dissipation range increases with concen- 
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FIG. 2: (color online) Concentration effects on the energy trans- 
fer rate et and energy dissipation rate ed at fixed Reynolds num- 
ber (R x = 350; Wi = 6.0). While e D (O) decreases smoothly 
as the polymer concentration increases, the energy transfer rate er 
(□) changes only when the concentration is above 5 ppm. u' 3 (A) is 
shown as a surrogate for the energy injection rate ei . The inset shows 
the change of u with concentration. The slow decrease of u is most 
likely due to the reduction of the forcing efficiency in the boundary 
layer of the propeller by drag reduction. 



tration. At higher concentrations ((f) > 7 ppm), the appar- 
ent inertial-range plateau is suppressed, indicating that the 
energy transfer rate has changed. We also measured the ef- 
fect on the energy dissipation rate eo by plotting in Fig. |TJb) 
£>Aw( r ) compensated by the dissipation-range scaling pre- 
diction (i.e., [2^D/vjv(r)/15r 2 ]). Due to the finite spatial 
resolution of our measurement system, the very small scales 
(r <C rj) are only partially resolved. This effect is more 
pronounced for pure water case since for polymer solutions 
the small scales increases. Within experimental uncertainty, 
for water, ez?(0) = er(0) as indicated by the dashed line in 
Fig. [Tfb). Nevertheless, we observe that the effect of poly- 
mers on dissipation scales is smooth, without any sign of a 
transition. 

To quantify the effect on the energy cascade, at each poly- 
mer concentration we measured er using Eq.[T]and tu using 
Eq.|2] as shown in Fig. [2] It is clear that et remains approx- 
imately unchanged for small concentrations, but drops sud- 
denly when the concentration is above 5 ppm. The energy 
dissipation rate en, on the other hand, decreases smoothly 
with concentration. We cannot measure the energy injection 
rate ej directly in our current apparatus. Since ej ~ u' 3 /L, 
however, measurements of the root-mean-square turbulent ve- 
locity can serve as a surrogate, provided that the integral scale 
L remains constant. The measured vf decreases slowly with 
concentration, as shown in Fig. [2] This nearly linear decay 
cannot account for the change of ey. 

The physical basis for the observed transition at a concen- 
tration of roughly 5 ppm remains unclear. It is possible that it 
may result from polymer-polymer interactions that occur only 



above the so-called "overlap" concentration. For the poly- 
mer used in our experiments, the overlap concentration based 
on the maximum extension length is ~ ppm, while it 

is approximately 200 ppm if based on the radius of gyration. 
Neither of these estimates coincides with the critical concen- 
tration observed in our experiments. 

In his pioneering work in the 1970s, Lumley suggested that 
the length scale of any polymer effect should be determined 
purely by the time scale on which the polymer recoils jjj. 
This (inertial-range) scale r* — (ctt^) 1 / 2 is the scale at 
which the local Weissenberg number is unity, and at which 
we expect the polymers to begin to be stretched by the flow. 
In this theory, there is no concentration effect, and so it cannot 
fully explain our experimental observations. It may, however, 
apply below the critical concentration. Since the scale r* in- 
creases with Reynolds number, we tested Lumley's hypothe- 
sis by varying the Reynolds number while keeping the con- 
centration fixed at 5ppm, just below the concentration where 
we observe changes in the energy transfer rate. The change of 
Dnn(i") with Rx is shown in Fig. a). As shown in Fig. Ob), 
however, scaling by r* does not collapse our data. We find in- 
stead, as illustrated in Fig. [3c), that the naive normalization 
of r by the Kolmogorov length scale rj (measured from pure 
water at the same Reynolds number) collapses all the curves. 
We note that rj decreases with Reynolds number in our ex- 
periments, while the prediction of r* increases. The behavior 
we observe is therefore qualitatively different from drag re- 
duction at a boundary, where Lumley's theory does appear to 
apply M. 

In contrast to Lumley's argument, Tabor and de Gennes 
suggested qualitatively that while the polymers are affected 
by the flow at r* , it is only at a smaller scale r** that the flow 
is affected by the polymers 112.71 12811 . In their framework, r** 
is determined by balancing the turbulent kinetic energy at a 
given scale with the elastic energy in the polymer phase. Each 
polymer molecule can store some elastic energy; increasing 
the number of polymers therefore increases the energy in the 
polymer phase. The Tabor-de Gennes picture thus allows for 
a concentration dependence of the polymer effect. There are, 
however, many undetermined parameters in their qualitative 
theory, and further development is needed to make a quantita- 
tive experimental test. 

In summary, we investigated the effect of very small con- 
centrations of long-chain polymers on the dynamics of tur- 
bulence in the bulk of the flow. We observed a strong con- 
centration dependence: the energy dissipation rate appears to 
change for all polymer concentrations, while the inertial-range 
energy transfer rate only changes above a critical concentra- 
tion. At a fixed concentration below the critical concentration, 
we observed that the scales at which the polymers affect the 
energy cascade decrease with increasing Reynolds number, in 
a manner similar to the dependence of the Kolmogorov scale 
rj on Reynolds number. Both the concentration effect and the 
Reynolds number effect can not be explained by the theory 
(using Lumley's time criterion) that has been used to explain 
the drag reduction phenomenon. Our results suggest several 
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FIG. 3: (color online) Reynolds number effects at fixed concentration. Compensated Eulerian structure functions are shown for four Reynolds 
numbers in a 5 ppm polymer solution, (a) The length scale of the polymer effects changes with Reynolds number, (b) Scaling by the Lumley 
scale r* does not collapse the data for the different Reynolds numbers, (c) The small-scale data collapse when r is scaled by rj, the Kolmogorov 
length scale determined from water data. 



challenges for future research. The qualitative change in the 
polymer effect above the critical concentration must be ex- 
plained, and the exact ways in which the polymers change the 
three energy rates in turbulence must be clarified. And finally, 
if Lumley's time criterion does not hold in the bulk, a physi- 
cal mechanism by which the polymers can affect scales much 
larger than their size must be identified. 
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